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ABSTRACT 

The unsteady response of two oceanic density fronts to 
local atmospheric forcing, using combinations of wind stress 
and surface heat flux, is investigated with an embedded 
mixed layer-general circulation model. The adjustment of 
the frontal structure is de pendent upon the wind stress 
direction and whether there is surface heating or cooling. 
In cases of an applied wind stress alone where denser water 
is transported toward less dense water, the frontal struc- 
ture diffuses, the mixed layer depth deepens, and cross- 
frontal mixing occurs. In cases where less dense water is 
transported toward denser water, the frontal structure is 
preserved, mixed layer depth is preserved and cross- frontal 
mixing is minimized. The addition of surface heating shal- 
lows the mixed layer and inhibits vertical mixing. Inertial 
oscillations are observed in the across- front velocity 
field. 
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I. INTRODUCTION 



A. EACKGRCUNE AND PURPOSE 

Curing recent years, there have been a number of scien- 
tific studies that describe the nature and existence of 
oceanic fronts. Cromwell and Reid (1956) identified a frcnt 
by an abrupt horizontal density change at the sea surface. 
Roden (1976) seated that a front should be considered to 
exist wherever any oceanic state variable (density, tempera- 
ture cr salinity) reaches a relative maximum value in its 
horizcntal gradient. 

The motivation for studying fronts and ocean prediction 
has been stressed by Roden (1976), Elsberry and Garwood 
(1979), and Niiler (1 982) . Elsberry and Garwood (1979) spe- 
cifically pcirt out antisubmarine warfare applications, 
fisheries management and a resource for climate research 
data, as viable and pertinent reasons for further study, 
.lany important physical processes which are accompanied by 
attendant air-sea interactions occur in frontal regions. 
Latent and sensible heat transfer, evaporative and precipi- 
tative processes, barcclinic currents, salt fluxes, mass and 
momentum transport, inert ial-int ernal wave activity and mix- 
ing all can occur. Oceanic frontal areas are areas of sound 
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speed changes and often are areas of increased biological 
activity. 

This thesis examines the time-dependent response of a 
numerically- simulated oceanic front to local atmospheric 
forcing. The model used in this work incorporates the forc- 
ings of wind stress and surface hearing with a 20-level 
primitive equation ocean circulation model with embedded 
mixed layer dynamics. 

B. OBSERVATIONAL STUDIES 

Cromwell and Reid (1956) described a simple case of an 
upper layer of warm, less dense water overlying a layer of 
cool, denser water in the ocean. The front in this case is 
a density front formed by the horizontal temperature gradi- 
ent. Observations shewed surface convergence along the 
front, sinking motion in the immediate vicinity of the front 
and divergence at lower depths. They observed that the wind 
affected the front by mixing the upper layer, reducing the 
horizontal surface temperature gradient, and enhancing the 
temperature jump at the base of the mixed layer. 

Voorhis and Hersey (1564) observed a wintertime thermal 
front in the Sargasso Sea. They found an average tempera- 
ture gradient of 1 deg C/10 km, and a maximum in one 
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instance cf 1.5 deg C/5 km. This was accompanied by a 60 
cm/s jet flowing to the east alcr.g the edge of the front. 
The depth of the front corresponded to the mixed layer 
depth. The water on either side of the front was found to 
be well-mixed. 

Katz ( 1969) continued the work of Voorhis and Hersey and 
derived temperature-salinity correlations for a front 
observed in the Sargasso Sea in May. He found the front to 
be a separation cf two distinct water masses, with a tran- 
sition zcne cf only a few meters in the vertical. He 

deduced that water transport along the front maintained the 
sharp transition. Although not observed, Katz presumed mix- 
ing must occur at the interface, but the along-front current 
prevented any build-up of well-mixed waters in that area. 
Ha also hypothesized that a vertical component to the cur- 
rent must exist at the interface. Thus, the interface would 
be preserved by both horizontal and vertical motions, and 
the vertical velocity would account for the observed accumu- 
lation of well-mixed waters below the front. 

Eang ( 1973) studied a front at the southern end of the 
Benquela Current. This was an intense front with a horizon- 
tal temperature gradient cf 1.6 deg C/km and a temperature 
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change of 8-10 deg C across it. There was also local 
upwelling occurring due tc an offshore transport. As might 
be expected from Ekman theory. Bang found the front to 
intensify under southerly winds (Ekman transport being to 
the left of the wind in the Southern Hemisphere) . Under a 
northwest gale, the front rapidly weakened. It reintensi- 
fied as the wind steadied from the south. 

Ecden has made significant and extensive contributions 
to the literature on oceanic fronts, especially on the 
large-scale fronts in the Pacific Ocean. In a series of 
papers, he has observed and examined the Pacific subarctic 
front (Roden 1975, 1977) , the subtropical front (Roden 1974, 
1975, 1980), the doldrura front (Roden 1974, 1975), and the 

subarctic-subtropical transition zone (Roden 1971, 1972). 

In the 1971 paper, Rcden concentrated on the subarctic-sub- 
tropical transition zone and found that the wind stress 
played an important role in determining the location of heat 
and salt flux divergence zones. Roden (1972) focused on the 
temperature and salinity fronts located at the boundaries of 
the transition zone; the subarctic front to the north and 
the subtropical front (which is also a density front) to the 
south. The origins of bcth fronts were found to be related 
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to the wind stress, geostrophic flow fields and heat and 
salt flux divergence. The doldrum front was included in a 
study of the subtropical front (Roden 1974) . Further work 
on all three fronts was conducted by Roden (1975) with ref- 
erence tc the wind and to the energy flux fields. He 
observed that the depth of the mixed layer changed dramati- 
cally across the subarctic front in winter and spring; from 
100 m on the north side tc 300 m on the south side. The 

subtropical front is highly dependent on the wind and energy 
flux fields at the sea surface, as the front occurs in an 
area of Ekman transport convergence. The doldrum front is a 
front due tc a salinity gradient only, with a baroclinic 
east ward- flowing current at the surface and a faster west- 
ward-flowing undercurrent. 

For the subarctic front, Roden (1977) found atmospheric 
forcing (that is, wind stress, radiative heat flux, precipi- 
tation and turbulent energy fluxes) to be the primary cause 
of frontal movement, intensification and decay. Roden has 
shown repeatedly that the roles of wind stress, Ekman trans- 
port and heating are inextricable from the dynamics of frcn- 
togenesis. Eoden (1980) again examined the subtropical 
frontal zone (within which are located a number of 
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individual fronts) and showed that Ekman transport led tc a 
concentration of the temperature and salinity gradients in 
the upper ocean. Thus, fron togenesis was due to horizontal 
convergence and confluence of the flow field. 

While Eoden's work has dealt primarily with frontogene- 
sis ard in explaining the dynamics of several identified 
North Pacific fronts, his findings are pertinent to this 
thesis, i.e., local atmospheric forcing plays a major role 
in the dynamics of oceanic fronts. 

C. MODEL STOGIES 

Garvine has concentrated on a modelling, approach to 
fronts. In one of his earlier papers ( 1974), he investi- 
gated the dynamics of small-scale fronts, though without 
regard to atmospheric forcing. Such fronts propagate into 
the ambient water via the unbalanced horizontal pressure 
gradient. Because of minimal cross-frontal mixing, he found 
strong convergence and sinking at the front. He later 
expanded this model to include the effects of wind stress 
(Garvine 1979a, 1979b) . However, because of the need to 

maintain steady state in the model, he prescribed wind 
stress to be of secondary importance relative to the hori- 
zontal pressure gradient. Buoyancy and energy budgets were 
added later tc the model (Garvine 1980). 
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Kao (1980) investigated the Gulf Stream as a large - scale 
density front in the upper ocean. The heart of his work was 
performing a scale analysis, which revealed three length 
scales: an inertial or deformation length scale, a buoyancy 

length scale, and a diffusive length scale. Kao’s conclu- 
sion was that the front was maintained by the cross-Gulf 
Stream circulation, but his model included no wind stress or 
heating. 

Cushman-Rcisin' s (198 1) model was based on Roden’s work 
on large-scale frontogenesis, and it does have some similar- 
ities to the model used in this thesis because it includes 
mixed layer dynamics. An important difference is in the 
scales. Cushman-Rois in investigates large-scale and long- 
term spin-up cf a front in response to a wind stress curl. 
The model in this thesis is used to investigate local 
respcnse of a pre-existing front to local winds and heating. 
Surface wind stress is a primary means of forcing in both 
models, but Cushman-Roisin applies a strong wind stress curl 
in the form cf a wind field that changes direction at the 
latitude cf the zonally-oriented front. Me, on the other 
hand, spin up a horizontally uniform wind stress over six 
hours, after which it maintains a constant value. The 
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results obtained by C ushman-Rois in depend much upon the curl 
of the wind stress. Cushman-R oisin finds the mixed layer 
depth is at a minimum in the middle of the front, while 
Ekman downwelling is at a maximum. This is probably due to 
a lack of any wind stress at the front itself, although the 
curl is large . 

Camerlengc (1982) studied the large-scale response of 
upper ocean fronts, specifically the Pacific subarctic 
front. His model has a variable-spaced horizontal grid with 
highest resolution in the frontal area. Vertical mixing was 
not included in the model. Assuming a pre-existing front, 
Camerlengc applied different wind stresses, and found that 
except in the case of the passage of a strong cyclone, the 
effects of wind forcing were limited to the upper 150 m. 
One case of a uniform wind stress did show Ekmar. transport 
and a horizontal displacement of the front. 
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II . DESCRIPT ION OF THE PRO B LEM 

A. GENERAL 

The response of a pre-existing oceanic front to local 
atmospheric forcing (wind stress and surface heating) over a 
time period of two days is examined. An embedded mixed lay- 
er-ocean circulation model is used, and two different frcn- 
tal structures are treated. The first numerically-inserted 
front, hereafter referred to as Front 1, is somewhat of an 
artifact in that the isotherms become horizontal at a rela- 
tively shallow depth (see Figs. 1,2, and 3 for the initial 
conditions) , so that Front 1 is manifested only within the 
upper 60 m. The temperature field is vertically stratified 
belcw approximately 70 m. The front is about 10 km wide. 
The maximum horizontal temperature gradient at the surface 
is 1.5 deg C/1.6 km. The mixed layer depth is a numerical- 
ly-inserted approximation to what the actual mixed layer 
depth would be. A geostrcphic along-front velocity is cal- 
culated from the initial temperature field. It has a maxi- 
mum speed cf 32 cm/s with its axis located at the surface on 
the warm side cf the front. The across-front geosrrophic 
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velocity is zero. Front 1 is a viable representation of a 
smaller mescscale upper layer front, and the model results 
will be important in this context. The second type of 
front, hereafter referred to as Front 2, is a simulation of 
a more typical front, similar to those reported by Roden 
(1980). The horizontal temperature gradient at the surface 
is 1.5 deg C/3.2 km. The front is not only evident at the 
surface, but also its effect on the temperature field 
extends throughout the thermocline of the model (see Figs. 
4, 5, and 6 for the initial conditions) . The mixed layer 
depth again is a numerically- ins erted approximation based on 
the shape of the temperature field. It is felt that adjust- 
ments which occur shortly after the start of integration of 
the model will compensate for any errors in the initial 
field. A gecstrophic alcng-front velocity is calculated 
from the initial temperature field. Its maximum value is 
100 cm/s and the axis lies at the surface on the warm side 
of the front. The acrcss-front gecstrophic velocity is 
zero. With both fronts oriented in a right-hand coordinate 
system, z increases upwards, y increases to the right, and x 
increases out of the page. 
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B. HYPOTHESES 

It is expected that the direction of the wind relative 
to the orientation of the front plays an important role in 
the ne&r-frcntal circulation and density structure- Ekman 
transport across the front should cause upweliing and down- 
welling features near the surface manifestation of the 
front. Johannessen , et al (1977) observed that the thermo- 
cline through the Maltese front shallowed in certain areas, 
which indicated upweliing, and also -chat the thermocline 
spatially deepened after isotherms surfaced. In Johannes- 
sen 's case, this was attributed to the presence of an eddy. 
If the water masses had been of similar types, it would have 
been indicative cf downwelling. 

With regard to one-dimensional mixing, Niiler's (1975) 
mixed layer model used an impulsive wind stress and no heat- 
ing, and showed that initially the mixed layer deepened rap- 
idly. Strong inertial motions can be generated which can 
lead to a large increase in available turbulent energy and 
thus to a rapid (within the first half pendulum-day) deepen- 
ing cf the mixed layer. Mellor and Durbin (1975) also used 
an impulsive wind stress of 2 cm /s (approximately 2 
dynes/cm^ ) and no heating. They showed that the mixed layer 
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T at hour 0 




Fiqure 1. Front 1 Initial Temperature. Only the middle 
section of the grid and the upper 175 m are 
shewn. Contour interval is 0.5 deg C. 
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H at hour 0 




Figure 2. Frcnt 1 Initial Mixed Layer Depth. 
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U at hour 0 




Fiqure 3. Front 1 Initial U-Geo strophic Velocity. Contour 
interval is 8 cm/s. 
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Figure 4. Front 2 Initial Temperature. Only the middle 
section of the grid and the upper 225 m are 
shewn. Contour interval is 0.5 deg C. 
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H at hour 0 




Figure 5. Front 2 Initial Mixed Layer Depth. 



26 



Depth (m) 



U at hour 0 




Figure 6. 



Front 2 Initial (J-Geo strophic. Velocity . Contour 
interval is 20 cai/s. Dashed lines mdrcate 
negative values. 
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deepened and that the temperarure jump also increased 
slightly. In a two-dimensional model, De Szoeke (1980) 
examined the effects of wind stress only in his version of 
the Niiler (1975) model and found that in the case of 
constant wind stress, the mixed layer deepened rapidly in 
the first half pendulum-day (t=Pi/f), and then slowed. An 
additional effect of the wind stress is the possible augmen- 
tation of a guasi-ge cstrcFhic frontal shear (in the aiong- 
front directicn) which could increase vertical mixing as 
prescribed by the initial value of the Gradient Richardson 
Number: 



_ s h. 

% bz 







2 




dz 





A surface heating-ccoling cycle can lead no stratifica- 
tion under light to moderate winds which will tend to modu- 
late any nixing tendencies brought about by wind-generated 
turbulence. Niiler (1975) showed that surface heating halts 
the deepening cf the mixed layer. We expect the same result 
in our model. Mellcr and Durbin (1975) investigated a case 
of a combined constant wind stress (2 dynes/cm^) and a 
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constant heating/cooling value (H = +.01 cm-K/s) . For the 
heating case they found a shallowing and leveling off of the 
depth of the mixed layer and an ever-increasing temperature 
jump value at the base of the mixed layer. The same heating 
condition with a lower wind stress (1 dyne/cm^) resulted in 
a shallower mixed layer depth and the establishment of a 
double thermocline, a remnant of the initial mixed layer and 
the formation of a new shallower layer due to the influence 
of heating. 

It is expected that inertial currents across the front 
may play an important role in the across-f rontal transport 
of heat f mass, and momentum. Phasing of inertial wind- 
driven currents in conjunction with vertical mixing may 
result in lateral entrainment or ’’capture 1 ’ of buoyant (warm) 
water by the adjacent cooler water during horizontal oscil- 
lations across the front plane. 



C. METHOD 

To test the above hypotheses and to establish the basic 
dependence of frontal response to local atmospheric forcing, 
four combinations of surface boundary conditions are inves- 
tigated. These are listed in Table I. 



29 



TABLE I 



Atmospheric Forcing used in Model Suns 



x (dyne s /c m ) 



(dy nes/ cm ) 



2 (d^2 C-cm/s) 



Case I +1.0 
Case II -1.0 
Case III +1.0 
Case IV -1.0 



0 

0 

0 

0 



0 

0 

-.004 

-.004 



In the cases above, T x is the along-front (x-direction) com- 
ponent of wind stress, is the cross-front (y-directict) 

component cf wind stress, and heat flux Q is applied to the 
ocean uniformly at the surface (heat flux is positive 
upwards, sc the negative value indicates warming of the 
ocean from the air) . All four cases are examined for both 
Front 1 and Front 2. 

Cases I and II examine the effects of wind stress paral- 
lel to the front only, which will eventually generate a 
steady Ekman transport perpendicular to the front. The wind 
stress of 1 dyne/cm^ corresponds to a wind speed cf 7-8 m/s. 

In all cases the wind stress is spun-up from 0 at t=0 to its 

2 

maximum absolute value (1.0 dynes/cm ) by hour 6. From hour 
6 to hour 48 it remains constant. Such a short spin-up time 
is net unreasonable, as winds can change on a shorter time 
scale, for example, when an atmospheric cold from passes. 
Cases III and IV impose a constant surface heating on each 
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of the two preceding cases. As indicated by the earlier 
one-dimensional mixed layer models, we expect that the hear- 
ing value will net be overwhelmed by the wind stress, and 
that the effects of both dynamic and thermodynamic forcing 
will be evident. 

The initial conditions for Front 1 are shown in Figs. 1, 
2, and 3. Warmer, less dense water resides on the left, and 
cooler, denser water resides on the right. The initial con- 
ditions for the velocity are determined by the associated 
initial density field, which is balanced by a geostrophic 
along-front current directed out of the figure (u positive) 
(see Fig. 3). The mixed layer depth varies from a uniform 
depth in the far field to a shallowing through the front. 
The model is run for a 48 hour time period for each of the 
specified atmospheric forcings shown in Table I. 

In Case I we expect tc see water from the cold (right) 
side transported to the warm (left) side, upwelling on the 
cold side at the front and dcwnwelling due to convective and 
turbulent mixing cn the warm side. The mixed layer depth 
should deepen, and the front should diffuse. In Casa II we 
expect the warm water to be transported on the surface to 
the ccld side, but convective mixing to be absent due to the 
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stable condition of the less dense water overriding the den- 
ser water. Upwelling should occur on the warm side at the 
front, and only minor downwelling should occur on the cold 
side due to the effect cf wind stirring alone. In both 
Cases III and IV the heating should have a shallowing influ- 
ence on the mixed layer depth tending to offset deepening 
and ether effects caused by turbulent mixing. 

The initial conditions for Front 2 are shewn in Figs. 4, 
5, and 6. In this configuration, the warmer, less dense 
water resides on the right and the cooler, denser water 
resides on the left. The accompanying along-front gecs- 
trcphic current is directed into the page (u negative) with 
its maximum value at the surface. The mixed layer depth 
follows the ccntcurs of the upper thermociine, shallowing 
across the front. The model is initiated and integrated in 
time as it was for Front 1. 

In Case I (wind stress positive, directed out of the 
figure) , we expect the Ekman transport to be toward the 
left. The results for this case should be analagous to the 
expected results for Case II applied to Front 1 as previ- 
ously noted. Case II should be compared to Case I for Front 
1 and Cases III and IV should be compared to Cases IV and 
III respectively cf Front 1. 
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The inertial period for our model is 





20 hours 



at the latitude cf 36 deg N. Considering that in this model 
wind stress dees not attain its steady maximum value until 6 
hours have elapsed, we should find a maximum in inertial 
motion somewhat later than half an inertial period after 
time zero, at perhaps t= 1 5 hours or sc. 



D. DESCRIPTION CF THE MODEL 

The primary documentation for the model is in Adamec, a: 
al (1981). In that paper, the model is cast in radial coor- 
dinates to investigate the oceanic response to a stationary 
ra dially-symmetric hurricane. In our case, the model is in 
right-hand Cartesian coordinates and it is two-dimensional 
(y and z) . It is believed that this model is particularly 
well-suited tc the study cf upper ocean fronts because it 
incorporates the Garwood (1977) bulk mixed layer model into 
the Haney (1980) primitive equation model. The ocean is 
assumed hydrostatic and incompressible, where density is a 
function of temperature alone. The Coriolis parameter is 
constant and there are no fluxes of mass, momentum or heat 
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normal to any boundary except at the surface, where atmos- 
pheric forcing is applied. The horizontal grid spacing is 
125 m and the vertical grid spacing is variable over 20 lay- 
ers from 6 m at the top layer to 64 m at the bottom layer. 
The dcmain size in y and z is 100 km by 500 m, respectively. 
With the front pre-existing as an initial condition, the 
initial u-ccmponent of velocity is in geostrophic balance, 
and the initial v-component is initially set to zero. 

Drawing upcn Adamec, et al (1981), highlights of the 
model are outlined below. Table II lists the pertinent 
model variables and constants used in the governing 
equations . 



The governing equations in their two-dimensional form are as 



follows: 
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TABLE II 

Model Constants and Variables 



<* 

at 


coefficient of thermal expansion 
horizontal eddy conductivity 


2. OE-5 


/degC 

p 


Am 


coefficient for heat 
horizontal eddy diffusion 


2.5E+5 


CII7S 

? 


JJU 

Km 


coefficient for momentum 
vertical eddy conductivity 


2. 5E+5 


cm/s 


i 

Kw 


coefficient for heat 
verzical eddy diffusion 


0.5 cm 


2 /s 


M 


coefficient for momentum 


0.5 cm 


2 /s 




latitude 


36 N 


f 


Ccriclis parameter 


8. 55 E- 


5/s 


g 


gravity 


981.0 


cm/s2 


?0 


reference density of seawater 


1. 0276 


g/cm^ 


T 0 


reference temperature 


5 degC 




At 


mcdel time step 


30 s 




J 


horizontal grid array size 


80 0 




K 


vertical and array size 


20 




Z=-D 


basin depth 


500 m 




y 

u 

V 

T 

h, H 


basin width 

x (alcn g-f r ont) component of 
velocity 

y (across-f ront) component of 
velocity 
temperarur e 
mixed layer depth 


100 km 





0 * - If - ?e (5) 
e - <? 0 c 1 - V > (6) 



The boundary conditions are: 



u = 0 





@ y = 0 , 100 km 



(7) 



w = 0 @ z = 0 , -D 



( 8 ) 
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-w'T’ 



Tx 

?o 



0 @ z 



0 




-u ' w ' 



-v'w' = 



-w 'T ' = -u ' w ' = -v'w' = 0 © z = -D 
An equation fcr the depth, h, of the well-mixed 

derived by integrating the continuity eguation and 
rhe rigid lid boundary condition (Eg. 8) : 



( 9 ) 

layer 



is 



applying 



bh. 

Ft 



+ w 



-h 



vv 



(io) 



The vertical turbulent fluxes are parameterized in two 
different manners depending upon where in the water column 
they are located. Below the mixed layer, the vertical 
fluxes are parameterized by eddy viscosity and conductivity 
coefficients : 



u ' w ' = 



V hi 
" dBz 



v'w' = 



-K & 
K M3z 



(ii) 






Above the mixed layer, turbulent mixing makes such a method 
which assumes and constant unrealistic and inaccurate. 
Thus the Garwood (1977) turbulence closure model is invoked 
using bulk turbulent kinetic energy equations. 
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The entrainment buoyancy flux is given by: 



e*g(w’T ' )_k = - <w 'w ') 2 (E)/h 



( 12 ) 



where < w * w •> aEd <E> are the vertical and total components 
respectively cf the turbulent kinetic energy. The mixed 
layer total turbulent kinetic energy equation and the verti- 
cal component of turbulent kinetic energy equation are: 



JL 



h<E>' 



= mu: 



*gh(w'T')_ h /(2 R i*) 




(13) 



2 f h(w 1 w 1 ) 
bt 2 



*gh((w’T’)_ h “ (WT*) 0 )/2 
+ ((E) - 3 (w 'w •) ) <2) 2 
- (<E>^ + fii)(E>/3 



(H) 



where u* = ( ) 2 



is the friction velocity and 



ai*= 



c< ghAT 

2 2 

Au + £>V 



is the bulk Richardson Number. 



The downward fluxes of momentum 
puts 4 as: )_ h . „ eiu 

-(v'w')_ h = » e Av 



- (u ' w ' ) and - (v ' w* ) 




are com- 



( 15 ) 
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where w 0 is the ent rainmert velocity: 



w 



e 



-(w'T')_ h 

ZTt 



( 16 ) 



The entrainment heat flux (w'T*) is the lower boundary 

-h 

condition for the mixed layer temperature profile, and a new 
mixed layer depth and new momentum and temperature profiles 
are calculated. The mixed layer depth is independent cf the 
model level depths. If the mixed layer shallows rather than 
deepens (as would occur if there were insufficient vertical 
turbulent energy to transport heat down to the depth cf the 
existing mixed layer), then the depth where (w'T 1 ) — > 0 
becomes the new mixed layer depth. During mixed layer shal- 
lowing, heat, momentum and potential energy are conserved. 
To prevent the levels beneath the mixed layer from becoming 
unstable in the numerical process, a dynamic stability con- 
dition is imposed such that the Gradient Richardson Number 
is always greater than 1/4. 

The embedding of the mixed layer model within the gen- 
eral circulation model requires communication between the 
two. The general circulation model is the dynamic portion 
which calculates changes of u, v, and T at each depth level 
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due to advective and diffusive processes 



These values 



are 



then acted upon by the m 
culatss the changes due 
mixing, and calculates 
values are then used by 



ixed layer model, which in turn 
tc surface fluxes and er.trai 
the new mixed layer depth, 
the general circulation model. 



cal- 

nment 

These 
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III. MODEL R3S0LTS FOR FRO NT 1_ CASS I 
A. ANALYSIS AND DISCUSSION 

Although the front is initially in geostrophic balance, 
the imposition of atmospheric forcing causes the front to 
undergo a variety of adjustments. In Case I, a positive 
wind stress cf 1 dyna/cm^ is applied. The front and the 
u-component field undergo diffusion, thereby giving rise to 
an ageostrophic u-component, which in turn creates a v-ccm- 
ponent. As diffusion occurs through mixing, and as the sur- 
face expression of the front is advected to the left in 
response to the wind stress, the along-front jet decreases 
in magnitude as well as in horizontal extent. A small 
countercurrent is present at hour 24 due to a reversal in 
the slope of the isotherms in the vicinity of the front. 
This effect disappears by hour 36, but then reappears at 
hour 48. 

The base cf the front on the right side deepens (as evi- 
dent by the isotherms becoming vertical to greater depths) 
from about 30 a initially to about 50 m as a result of the 
effect of turbulent mixing and the creation of ageostrophic 



velocities while undergoing dynamic adjustment (see Fig. 7). 
An ageostrcphic u-component is evident in Fig. 8 in the area 
of y = 53 Jem to 61 km. The mixed layer depth increases as the 
upward slcping thermccline erodes within the front because 
water is transported across the front, right to left. The 
mixing of the coder water with the warmer water causes the 
mixed layer depth to deepen on the left (warm) side. The 
axis of maximum v-values lies just above the bottom of the 
mixed layer near y=5 3 km. If the v-profile is examined 
(Fig. 9) , it is seen that dv/dy > 0 on the right (cold) side 
of the frontal axis, and dv/dy < 0 on the left side. The 
gradients of v are much stronger on the left side of the 
front than they are on the right. The maximum v-velocity is 
-8 tc -10 cm/s. From continuity, aw/dz < 0 (downwelling) on 
the right side of the axis and dw/dz > 0 (upweiling) in a 
thin layer on the left side that closely follows the mixed 
layer depth. It is noted also that the mixed layer depth 
has deepened at the right side of the front. This phenom- 
enon is present throughout the model integration. Though 
such a condition was not anticipated prior to the running of 
the model (we expected an upweiling in the mixed layer on 
the right side as a result of the mass transport over the 
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front towards 



rhe left) , this depression in deptn is 
believed tc be physically realistic. It is suspected that 
this increase in the mixed layer, depth occurs in response to 
no n- stationary adjustment to the atmospheric forcing, as it 
also occurs in Case II where the oppos it sly- directed wind 
stress is applied. Because the mixed layer depth responds 
to the vertical integral of transport, the fact that the 
mixed layer is shallower within the front than in the far 
field blocks the flow from penetrating to depth. As the 

water encounters this "obstacle", it "backs up" on the right 
side and thus forces a depression in the mixed layer. This 
explanation must be tempered by the fact that these figures 
are instantaneous snapshots in time, and can be deceiving in 
attempting tc describe processes that have been evolving 
over a length of time. Another possibility is that some 
type of internal standing wave has been excited which is 
causing fluctuations in the mixed layer depth. 

The influence of the ageostrophic u-veiocity, the core 
of the v-velocity lying along the well-mixed layer and the 
depression of the mixed layer on the right side of the front 
is felt at depths greater than 100 m. Upwelling and down- 
welling are evident at these depths (Fig. 7) beneath the 
right edge of the front. 
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Depth (m) 



T at hour 12 




Fiaure 7. Front 1 Case I 12-Hour Temperature. Contour 
interval is 0.2 deg C. 
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U at hour 12 




Figure 8. Frcr.t 1 Case I 12-Hour U-Velocity. Contour 
interval is 8 cm/s. 
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V at hour 12 




Figure 9. Front 1 Case I 12-Hour V-Velocity. _ Contour 
interval is 2 cra/s . Dashed lines rndicare 
negative values. 



45 



Interesting patte 
24 and 48. At 24 hou 
contrary to the expe 
surface layer, and i 
created by the inerti 
48 hours, two cores 
symmetrically centsre 
either side cf the 



fr on 


t, b 


ut form 


outs 


far 


field 


where 


the m 




The 


across- 


front 


with 


the 


d irecti 


on o 


ex ce 


pt ho 


ur 24. 


As 


the 


base 


of the 


mixed 


is n 


cn-ze 


rc ben 


eath 


fron 


t as 


well . 


Thu 


does 


occu 


r in th 


is ca 


vert 


ical 


mixing 


cf he 


B. 


SUMMA 


BS -- c 


ASE I 




Under 


the in 


fluen 


ing. 


the 


front 


ha s 



direction of the surf 



ms in the 
rs a sraa 11 
cted Ekman 
s probably 
al motion, 
of a posit 
d about the 
front. Th 
ide the fro 
ixed layer 
velocity 
f the net S 
noted, this 
layer rath 
the mixed 
s, it appe 
se, and i 
at in the t 

ce of posit 
diffused an 
ace Ekman t 



v profiles evolve at hours 
positive v exists. This is 
transport direction in the 
the result of an oscillation 
Furthermore, at both 24 and 
ive v-component are present 
base of the mixed layer on 
ese are not present in the 
ntal area and extend to the 
is relatively stable, 
component is in consonance 
kman transport at all times 
velocity is at a maximum at 
er than at the surface. It 
layer in the middle of the 
ars that interfacial mixing 
t may be contributing to the 
hermocline at this location. 

ive wind stress and no heat- 
d has been advected in the 
ransport. Because the front 



46 



is unsteady under an applied wind stress, ageostrcphic 
velocities are created. Fart of the response is frictional 
due to the turbulent boundary layer or mixed layer pro- 
cesses, and part is inertial due to the applied wind stress 
boundary ccnditicn. A noticeable deepening of the mixed 
layer occurs at the front on the right (cold) side, and the 
horizontal variability in mixed layer depth and temperature 
are reduced significantly under the influence of both verti- 
cal and crcss-frontal mixing. 
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IV. 



KCDEL SES CLTS FO R F RONT J. CASE II 
A. ANALYSIS AND CISCUSSICN 

In this case, a negative wind stress (directed into the 
page) is applied. Wind stirring must, overcome static sta- 
bility as the warmer water will override the cooler water. 
It will be seen that the wind stress cf 1 dyne/cm^ is suffi- 
cient to mix the warmer water which is transported toward 
the cooler water. There is no evidence of stratification of 
shallow waters. Recall from Case I that, as the denser 
wafer from the right side was being transported to the left, 
it had to overcome static stability in order to upwell and 
override the front, and that the mixed layer depth profile 
posed a vertical barrier to the transport. In this case, 
the base cf the mixed layer is at an incline xn the direc- 
tion of transport and static stability need not be overcome 
in the transport process. This leads to the expectation 
that cross-frontal mixing may be minimal in this case. 

Diffusion in the temperature field again is evident over 
the 48 hours, but this time it occurs primarily at the edges 
of the front. A strong horizontal temperature gradient of 
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0.6 deg C/1.2 km in the middle of the front is maintained up 
through 48 hours (see Fig. 10 as representative of this) 
which was not present in Case I. This lends credence to the 
hypothesis that cross-frontal mixing is minimal. 

Noticeable mixing occurs on the cold side of the front. 
The isotherms are vertical to greater deprhs than they were 
in Case I (see Fig. 10) . The turbulent mixing alone not 
only overcomes static stability but is also able to mix to a 
relatively greater depth on the right side of the front. 
The water than is being transported across the front is 
apparently immediately mixed by wind stirring as no plume of 
the warm, less dense water is present cn the right side. 

The alcng-front current again forms an ageostrophic com- 
ponent as it undergoes adjustment, and is especially evident 
on the right side of the front around y=55 km (Fig. 11) 
where it exists to a greater extent both horizontally and 
vertically than was seen in Case I. The core of the alcng- 
front jet dees decrease in speed and diffuse in time. By 
hour 36 the maximum value of u has been reduced to slightly 
more than 8 cm/s. The creation of this ageostrophic veloc- 
ity accompanied by a v-compcnent leads to strong downwelling 
around y=55 km underneath the front (Fig. 10) and remains 
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Depth (m) 



T at hour 36 




Fiqure 10. Front 1 Case II 36-Hour Temperature. Contour 
interval is 0.2 deg C. 
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over the 48 hour integration. A strong core of negative 
across-front velocity of ever -6 cm/s is created in this 
area (Fig. 12), and there is an attendant sharp upwelling 
in the isotherms at 75 m and wave-like fluctuations present 
to depths beyond 175 m. In the areas where dv/dy > 0 (dw/dz 
< 0) , downwelling is seen in the front, and where dv/dy < 0 
(dw/dz > 0) , upwelling is seen. In the surface layer, the 
v-component is positive and in the direction of Ekman trans- 
port. It is concentrated in a jet in the area of the tight- 
est horizontal temperature gradient over the shallow peak in 
the mixed layer depth. All of these tendencies and trends 
persist over the 48 hour integration. 

A deepening in the mixed layer occurs to the right of 
the front near y=55 km and progresses outward in time. The 
shallow peak in the mixed layer depth within the front is 
mostly preserved, though some erosion and mixing obviously 
has occurred (compare Figs. 13 and 14). This further sup- 
ports the hypothesis that minimal cross-frontal mixing 
occurs. Examination of the corresponding v-fields also show 
little crcss-frcntal mixing. 
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B. SUMMARY -- CASE II 

In this case of a forcing of a negative wind stress and 
no heating, the temperature structure of the front retains a 
strong horizontal gradient over the 48-hour integration. 
Diffusion of the front occurs, but it is concentrated at the 
boundaries of the front. Kind stirring overcomes the static 
stability of transported water once it crosses the front. 
The creation under adjustment of a large ageostrophic u-ve- 
locity and an accompanying v-velocity give rise to strong 
upwelling and dcwnwelling features within the front. The 
shallow peak cf the mixed layer depth within the front is 
preserved ever time. The fact that this peak is maintained, 
the horizontal temperature gradient retains a tight struc- 
ture, and the v-field is mainly constrained to the agecs- 
trophic areas and surface layer, suggests that interfacial 
or cress- frontal mixing is minimal. 
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U at hour 36 




Figure 11. Front 1 Case II 36-Hour O-Velocity. Contour 
interval is 4 cm/s. 
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V at hour 12 




Figure 12. Front 1 Case II 12-Kour V-Velocity. Contour 
interval 2 cm/s. 
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H at hour 24 




Figure 13. Front 1 Case II 24-Hour Mixed Layer Depth. 
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H at hour 48 




Figure 14. Front 1 Case II 48-Hcur Mixed Layer Depth. 
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V. 



MODEL RES UL1S FOR FRONT J. C ASE III 



A. ANALYSIS AND DISCUSSION 

A uniform surface heating is applied with the identical 
wind stress cf Case I. The magnitude of each forcing is 
such that influences from both are expected to be mani- 
fested. In examining Figs. 15 and 16, the obvious effect of 
heating is the adjustment cf the temperature structure on 
either side cf the frcnt. At first examination, it appears 
that the heating is having the opposite effect of what is 
expected - that is, it appears to cool the mixed layer 
rather than warm it. In reality, the heat is carried dcwn 
to depths cf 50 m and warms the water there. Cooling occurs 
between 32 and 50 m depth, which is readily seen in Fig. 17. 
Mixing may create a dynamic instability in the thermocline, 
and vertical diffusion of heat by Kq^d^T/dy^ may cause the 
intermediate water to be ceded and the deep mixed layer to 
be warmed. This basic structure remains over the 48-hcur 
integration. The frcnt undergoes diffusion and is advected 
to the left, as occurred in Case I. The front is manifested 
to a lesser depth in accordance with the heating process - 
downwelling and mixing have been moderated by the heating. 
Though seme wave-like motion is seen in the temperature 
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field beneath the mixed layer, the magnitude is not as great 
as in the case of wind stress alone. The mixed layer depth 
has shallowed tc about 43 m under heating alone in the areas 
outside the front (see Fig. 16). This shallowing occurs 
because the atmospheric forcing (heating and wind) is input 
as a body force ever the entire mixed layer. The heating 
has the effect of making the mixed layer depth more uniform. 
The marked deepening in the mixed layer depth that occurred 
in Cases I and II at the right (cold) side of the front is 
not present in this case. The deepening of the mixed layer 
affects the u-field. The agecstrophic u-compone.nt does not 
develop to as great an extent as it did in Casa I. The gen- 
eral along-frent current appears much as it did in Case I, 
except that its vertical extant is also constrained by the 
shallower mixed layer depth. In both cases, maximum speeds 
agree ever the 48-hour integration, and they are advacted 
with the front as it responds to the mass transport. 

The v-field follows the same trends as it did under wind 
stress alcne. At 12 hours, v is directed to the left as 
expected following the transport, and at 24 hours the v-ve- 
locity reverses and establishes the layers of concentrated 
positive v-values centered symmetrically about the mixed 
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T at hour 12 




Figure 15. Frcnt 1 Case III 12-Hcur Temperature. Contour 
interval is 0.2 deg C. 
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H at hour 12 




Figure 16. Front 1 Case III 12- Hour Mixed Layer Depth. 
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Figure 17. 



Front 1 Case III 
taken at y=35 km 



- T vs z Profile. Profile 
for t=0 and t = 1 2 hours. 
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layer depth. At 36 hours, v returns to a negative value ana 
at 48 hours, though still negative, v is much decreased in 
magnitude. Again, cells cf positive v-velocities are pres- 
ent at about the mixed layer depth. The across-front veloc- 
ity dees penetrate the front and the shallow peak of the 
mixed layer (see Fig. 18) which is indicative of cross-frcn- 
tal mixing. The shallow section of mixed layer itself is 
eroded and smoothes cut ever the 48-hour integration as it 
did in Case I (Fig. 19). 

B. SUMMABY -- CASE III 

The effects cf beth wind stress and heating are clear in 
this case. The wind stress induces a mass transport to the 
right, moving cooler, denser water into warmer, less dense 
water. This water not only apwells and moves over the 
front, but also moves across the frontal interface. At 24 
and 48 hours, the v-field appears to be influenced by 
inertial oscillations with a reversal in direction occurring 
in the surface layer at 24 hours and a near-reversal at 48 
hours. Cores cf positive velocities are established cen- 
tered about the mixed layer depth outside of the front. The 
front both diffuses and meves to the left in time. 
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Figure 18. Front 1 Case III 36-Hour V-Velocity. Contour 
interval is 2 cm/s. 
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Figure 19. Front 1 Case III 48- Hour Mixed Layer Depth. 
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Heating shallows the mixed layer depth as expected, and 
thereby restricts the vertical extent of the velocity 
fields. Heating also prevents the immediate deepening of 
the mixed layer on the right side of the front, which was 
also seen in tcth Cases I and II. The heat manifests itself 
at the initial mixed layer depth where the isotherms become 
more stratified. 
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VI. MODEL RESULTS FOR FRON T _1 CASE IV 
A. ANALYSIS AND DISCUSSION 

Case IV combines a negative wind stress and a uniform 
surface heating. The mixed layer shallows noticeably and 
heat is carried to the bottom of the mixed layer at depths 
of 55 to 70 m as in Case III (see Figs. 20 and 21). The 
heating dees have an interesting effect when compared with 
the wind- stress-only forcing of Case II. The sharp dip in 
the isotherms in the middle of the front is not evident here 
as it was in Case II. Horizontal diffusion of the tempera- 
ture structure occurs more rapidly and attains a gradient of 
1.6 deg C/5.6 km by hour 12 (see Fig. 20) . The along-frent 
velocity is limited in its vertical development on the right 
side of the front near y=57 km (Fig. 22) and follows both 
the horizontal and vertical diffusion tendencies of the den- 
sity front. The across-front velocity has a narrow, verti- 
cally-oriented core of negative values with a maximum of -6 
cm/s on the right side of the front between 25 and 50 m 
depth as in Case II. Since this maximum of v- velocity 
exists in an area where there is little u-velocity, its 
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existence cannot be attributed to a response to the agecs- 
trophic u-ccrapcnent. The position of this v-velocity corre- 
lates well with the upwelling/do wnwelling pattern evident in 
the isotherms in Fig. 20 , and also with wave patterns in the 
isotherms over the entire depth extent in the right half of 
Fig. 23. We suspect that such a large area of v-component 
created below the mixed layer is a response to the warming 
of the water at the initial mixed layer depth which has 
shoaled, re-oriented the isotherms and created vertical 
velocity. Furthermore, wind -generated mixing forces verti- 
cal motions in the water mass to a great depth, inducing a 
v-component there. 

In this case, the front diffuses horizontally more rap- 
idly than under wind stress forcing alone. The front 
spreads like a warm, less-dense plume overriding the cccl, 
denser water (compare Fig. 24 to Fig. 20) . The along-frcnt 
current shifts with the front. The mixed layer loses its 
shallow section within the front and smoothes out, which did 
not happen in Case II. The peak in the mixed layer was 
maintained much longer. With heating, there is less of a 
temperature difference and less of a shallowing from the 
f a r- field mixed layers to the shallow mid-frontal mixed 
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Figure 20. Front 1 Case IV 12-Hcur Temperature. Contour 
interval is 0.2 deg C. 
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Figure 21. Front 1 Case IV 12-Hour Mixed Layer Depth. 



69 



Depth (m) 



U at hour 12 




Figure 22. Front 1 Case 17 12-Hour 0-Velocity. Contour 
interval is 4 cm/s. 
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Figure 23. Front 1 Case IV 12-Hcur V-Velocity. Contour 
interval is 2 cm/s. 
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Figure 24. Front 1 Case IV 36-Hcur Temperature. Contour 
interval is 0.2 deg C. 
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B. SUHMABY — C A S3 IV 

Application cf a negative wind stress and surface heat- 
ing results in the front moving as a warm overriding plume 
to the right. The mixed layer initially shsliows due to the 
heating, and ever the 48 hour integration breaks down and 
smoothes cut in the horizontal. The along-front velocity is 
advected and diffused in the mixed layer. The pattern cf 
across-frcnt velocity shifts to the right in the surface 
layer, with its maximum value centered over the shallowest 
portion of the mixed layer. A v-ccmponent develops to quite 
a great depth on the right side of front due to motion gen- 
erated by heating and turbulent mixing. Cells of negative 
v- componen t are formed at hours 24 and 48 which are centered 
about the mixed layer depth. 
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VII- MODEL RESULTS F OR FRON T 2 CASE I 
A. ANALYSIS AND DISCUSSION 

The initial conditions for Front 2 are shown in Figs. 4, 
5 r and 6. The front extends throughout the deprh of the 
basin and is oriented such that warmer, less dense water is 
on the right and cooler, denser water is on the left. Many 
of the dynamical features of Fronz 1 are expected to be 
repeated fcr Front 2. In comparing Front 2 to Front 1, the 
horizontal axes are oppositely-directed to the frontal ori- 
entation. The ther mocline slopes upward in the positive 
y-directicn for Front 1, hut downward for Front 2. 

The wind stress cf Case I applied to Front 2 creates a 
transport directed to the left. Warmer, less dense water 
moves towards the cooler, denser water. In this case, dif- 
fusion in the temperature field is present for the duration 
of the integration. The horizontal temperature gradient 
remains large (1.0 deg C/1.2 km) in the middle of the front 
with the majority of the diffusion occurring at rhe right 
boundary. Mixing on the cool (left) side of rhe front 
extends frcm the surface to 135 m which is below the initial 
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depth of the mixed layer. Results of downwelling are seen 
at the left frontal boundary, near y=42 km. Apparently, 
there is no chance for a stable, overriding warm plume to 
develop before it is mixed into the ambient cool water peel. 
These phenomena in the temperature field were also seen in 
Case II fer Front 1. Much more dramatic oscillations occur 
in the isotherms, primarily on the right side between 53 km 
and 70 km, where the entire field is inclined. They are 
also visible or. the left side superimposed upon an inclina- 
tion of the entire field sloping downward to the right (Fig. 
25). The isotherms show a definite trend towards becoming 
more uniformly inclined beneath the mixed layer across the 
entire horizontal expanse, with the steeply-inclined field 
on the right becoming more horizontal. Downwelling causes 
the left side to slope downward to the right and upweliir.g 
lifts the isotherms on the right side. 

The mixed layer depth undergoes major changes over the 
entire domain. Cn the left side, it deepens to where the 
depth in the far field is approximately 125 m (Fig. 26) . 
Since this deepening occurs by 3 hours of integration, such 
a drastic deepening probably resulted from the initial mixed 
layer depth being inserted at too shallow a depth in this 
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Fiaure 25. Frcnt 2 Case I 36-Hour Temperature. Contour 
y interval is 0.5 deg C. 
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area. As in Front 1 Case II, a downward bulge in the mixed 
layer extending to 155 m occurs at the left frontal boundary 
at y=U3 km (Fig. 26) . In this case, the downward bulge dees 
not correspond to the side on which the waters would tend to 
"back up" prior to crossing over the front as in Front 1 
Case I. As strong mixing is implied by the vertical orien- 
tation of the isotherms at that location, the bulge is 
believed to he associated with a cooling of the waters in 
that area as a result of enhanced vertical mixing caused by 
the ageostrophic along-front velocity shear (Fig. 27) . 
Within the front, a shallow peak in the mixed layer depth is 
maintained at the left boundary, though over time it does 
deepen slightly and erode through mixing. There is much 
variability in the mixed layer depths on the right side of 
the front, but it appears that the general trend is for 
deepening to occur in the far field beyond y=55 km (Fig. 
26) . 

Along with the isotherms, the along-front u-velocity 
diffuses and decreases its maximum speed over time. As 
already seen, an ageostrophic component forms on the left 
frontal boundary and stimulates mixing there. The remnants 
of a tight horizontal velocity gradient are seen after HQ 
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Figure 26. Front 2 Case I 36-Hour Mixed Layer Depth. 
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Figure 27. Frcnt 2 Case I 12-Hour 0-Velocity. Contour 
interval is 1C cm/s. 
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hoars, though little cr nc advecticn of the along-front cur- 
rent occurs in response to the wind-generated transport. 
Wave-like oscillations appear to propagate along the deeper 
extremities cf the u-field, but no phase relation between 
isotachs or with ether fields is readily apparent. 

The acrcss-front v-velocity does correlate with the 
direction cf transport up through hour 9. An oscillatory 
pattern exists in the v field over the 48-hour integration. 
At hour 12, the entire v-field is positive, but this may be 
a response to the agecstrophic u-velccity, and it is strong 
enough to overcome the acrcss-frcnt transport velocity. At 
24 hours, v is in the direction cf transport and a maximum 
of -20 cm/s appears at the surface within the front over the 
shallowest section of mixed layer depths at y=43 km (Fig. 
28) , and remains through hour 36 . This is reassuring as it 
was expected that the inherent static stability of the water 
masses would minimize deep cross -frontal mixing, similar to 
what was seen in Front 1 Case II. The isotachs very defi- 
nitely terminate at the left boundary of the front, but they 
also exist beneath the mixed layer and do cross the front 
beneath the mixed layer. At hour 48, v is positive again, 
perhaps due to inertial motions. It is proposed (though it 
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is certainly not as clear as it was for Front 1) that 
cross-frontal nixing, though inhibited, occurs here tc a 
greater extent than it did for Front 1. The majority of the 
mass transport appears to exist in the mixed layer and 
passes over, rather than through, the shallow mid-frontal 
depth of the nixed layer. The reasons for the variations 
between the fronts may be due to the greater u-component 
(both in speed and in areal coverage) present with Front 2. 

B. SUMMARY — CASS I 

Several comp lexities have been incorporated into the 
modeling of Front 2. First, the isotherms slope downward to 
the right ever the entire right half of the basin. Hence, 
the front not cnly exists at the surface but is "connected” 
through the depth of the basin. This establishes an initial 
u-velccity which exists over the same area. Adjustment to 
the wind stress may not be confined solely to the upper few 
meters, as dynamic responses may be transmitted to depths 
beyond those directly affected in the surface layer. Also, 
the mixed layer depth is artificially inserted as an initial 
condition, and may net be the most perfect fit to the temp- 
erature field. Though this mixed layer configuration may 
exaggerate cr suppress actual features, it is nevertheless 
felt to be a reasonable approximation of an actual fronr. 
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Figure 28. Front 2 Case I 24-Hour V-Vslocity. Contour 
interval is 5 cm/s. 
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The responses are much more complicated in structure in 
Front 2 than they were in Front 1. Beyond 12 hours, the 
front shows little advective response to the mass transport. 
One reason this occurs is because the front exists over such 
a great depth. The original thermal structure and tempera- 
ture gradient are maintained in the middle of the front, 
with horizontal diffusion occurring primarily at the right 
boundary. Mixing in the surface layer occurs as the isot- 
herms become vertical to greater depths over time. At the 
left houncary of the front, downwelling is quite evident and 
the waters mix easily due to the creation of an ageostrophic 
u-component which enhances the wind stirring. Wave-like 
oscillations are evident in the temperature field beneath 
the mixed layer, and the entire deep temperature field 
becomes mere uniform in slope ever the 48-hour integration. 

The mixed layer depth must be considered in three sepa- 
rate regimes - left side, frontal, and right side. On the 
left, the mixed layer depth increases by almost 50 m at the 
frontal boundary, where mixing is enhanced due to the influ- 
ence of the acecstrophic current. The mixed layer bulges 
downward, then rises sharply to its shallow peak at the 
front. This shallow section remains throughout the 48 hour 
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integration though it dees deepen slightly due to mixing. 
On the right side, the mixed layer deepens and becomes more 
horizontal though there are strong fluctuations over inter- 
vals of 4 km and less. 

Two primary trends are observed in the along-front u-ve- 
locity: 1) the u-component diffuses in time, still retain- 
ing its basic shape but decreasing in speed and horizontal 
gradient; 2) an ageostrophic component is evident, prima- 
rily extending outwards from the left boundary of the front 
and which is believed to "drive" the downward bulge in the 
mixed layer there. 

The across-front velocity component is concentrated at 
the surface in the direction of transport between 0 and 9 
hours and at 24 and 36 hours. The front poses no barrier to 
the v-velccity at 12, 24, or 48 hours. The 36-hour fields 
are the only cr.es which have essentially no cross-frcnt 
velocity component other than in the surface layer. Because 
the initial horizontal temperature gradient is preserved in 
the middle of the front, and because the shallow portion of 
the mixed layer depth is still evident after 48 hours, sug- 
gest that little cross-frontal mixing occurs. 
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VIII. MODEL RE SULTS F OR FRONT 2 CASE II 
A. ANALYSIS AND DISCUSSION 

With the wind stress opposite from that of Case I, the 
net Ekman transport is left to right, and cooler, denser 
water will be transported toward warmer, less dense water. 
There are similarities here to the results of Front 1 Case 
I. The front bscomes much more diffuse at the surface and 
the horizontal temperature gradient is not preserved as it 
was in the preceding case (see Fig. 29). It has decreased 
to 1.5 deg C/6.4 km. The horizontal field on the left side 
reorients by inclining down to the right, and the inclined 
field on the right becomes more horizontal. The total 
effect is to make the temperature field more level beneath 
the mixed layer. Wave-like fluctuations are also seen in 
the isotherms (Fig. 29) . Vertical mixing can be readily 
seen to have occurred within the front by the steep slope of 
the isotherms. 

The mixed layer depth immediately adjusts in the left 
far field tc arcund 125 m as it did in Case I, but a much 
greater downward bulge in the mixed layer forms at the left 
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F-i aure 29. Front 2 Case II 36-Hour Temperature. Contour 
9 interval is 0.5 deg C. 
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boundary cf the front around y=4 2 kra and remains over the 48 
hour integration (Fig. 30). This result is consistent with 
the interpretation proposed for Front 1 Case I, where the 
shallower mixed layer within the front acted as a barrier to 
mass transport, causing the transported water to "back up". 
The mixed layer peak within the front deepens and erodes. 
The right side far field is again variable but shows the 
same tendencies to deepen. 

The agecstrophic u-velccity (not shown) is very similar 
to the isctach pattern which evolved for Case I. A core of 
maximum velocity is advected to the right which correlates 
with the advecticn of the front. The biggest difference is 
that the isotachs appear at a greater depth just inside the 
left frontal boundary. 

The expected direction of the across-front velocity is 
left to right, the direction of the net Ekraan transport. As 
in Case I, a wave oscillation is evident. At 12 and 48 
hours this is clearly so, with the velocity field existent 
along the entire section cf the front from the base of the 
front at around 150 m up to the surface (Fig. 31). The max- 
imum velocity occurs at a depth beneath the shallow mid- 
frontal section cf the mixed layer, which seems to 
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Figure 30. Front 2 Case II 48-Hour Mixed Layer Depth. 
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contribute to cross-frontal mixing, erosion of the mixed 
layer depth within the front, and weakening of the frontal 
horizontal temperature gradient. At 24 hours, the v-velcc- 
ity field is reversed, and negative velocities predominate 
in the frcrt and beneath the mixed layer. Inertial motions 
probably account for this reversal. 

B. SUMMARY -- CASE II 

Ihe temperature structure of the front shows vertical 
mixing and hcrizontal diffusion occurring over the 48-hcur 
integration. The horizontal temperature gradient is 
decreased due to the effect of convective mixing of the 
cooler water which is transported to the warmer water by 
wind stirring. The across-front velocity patterns penetrate 
the mixed layer within the front and establish cross-frontal 
mixing. The entire temperature field in the vicinity of the 
front undergoes an adjustment which inclines all the isot- 
herms downward tc the right in a rather uniform pattern. 

The mixed layer undergoes major adjustment over the 
48-hcur integration, but the primary trends are a deepening 
on bcth sides cf the front, a leveling off on the right side 
and a breakdown cf the shallow section within the front due 
to mixing and diffusion. By extrapolating the processes 
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Figure 31. Front 2 Case II 48-Hour V-Velocity. Contour 
interval is 5 cm/s. 
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occurring at 48 hours, it is presumed that the mixed layer 
depth would eventually even out with little or no remains of 
a frontal peak. As in other cases, a marked deepening in 
the mixed laysr occurs at the left boundary of the front. 
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IX. MODEL RESULTS FOR FR ONT 2 C ASS I II 
A. ANALYSIS AND DISCUSSION 

Constant surface heating is combined with a positive 
wind stress. As the net Ekman transport wiil be from right 
-co left, and with heating, we expect to see the formation of 
a warm plume cf water on the surface in the cool (left) 
ambient pool. The surface heating should enhance the inher- 
ent static stability and prevent the wind from completely 
mixing the transported warm waters. This effect does happen 
prior to the 12 hour point and becomes more pronounced over 
the entire integration (see Fig. 32). Downwelling and mix- 
ing cf the warmed surface waters with the cooler underlying 
waters are observed at the leading edge of the progressing 
front. Upwslling on the right side of the cooler water 
which replaces water lost due to transport is clearly seen 
by comparing successive temperature fields. The front shews 
a tendency tc almost split at 30 ra depth into a shallow sur- 
face front and a deeper front. Perhaps a different combina- 
tion of forcing would make this happen. 
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Figure 32. Front 2 Case III 36-Hour Temperature. Contour 
interval is 0.2 deg C. 
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The isotherms are much smoother than for the cases of 
wind stress alone. As in Case I, the originally- horizontal 
field on the left side becomes inclined, and the 
originally-inclined field on the right side becomes level. 
The effect is net as dramatic as before. 

The mixed layer quickly responds to the applied heating. 
By 12 hours the mixed layer has re-formed in the far field 
at 40 m. The shallow section within the front is advected 
with the plume of the front and moves into rhe cooler (left) 
side. The original shallow section of rhe mixed layer depth 
at the right boundary of the fronr has deepened no the mean 
far- field value prior no heur 12, "filling in" at the rear 
as the shallow frontal peak is advected. This shallow "bub- 
ble" progresses to the left with the front causing shallow- 
ing in advance of the front and deepening behind. The mixed 
layer responds with changes of up to 10 to 15 m on a time 
scale of something less than 12 hours (compare Fig. 33 tc 
Fig. 5). 

The along-frent velocity maintains a constant shape and 
configuration over the 48-hour integration. The maximum 
speed decreases slowly over time from 100 cm/s at the start 
to 6C cm/s by hour 36. The horizontal gradient on the right 
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Figure 33. Front 2 Case III 36-Hour Mixed Layer Depth. 
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side decreases rapidly but on the left side is preserved 
(Fig. 34) This is due to the advance of the plume and the 
resistance to mixing that the heating has established. 

At 24 and 36 hours, the across-front velocity is neg- 
ative, with the maximum of -10 cm/s occurring within the 
moving plume cf overriding warm water (Fig. 35) . A wave 
oscillation in the v field is evident, and has the same 
direction and time scale as these of Cases I and II. At 12 
hours, v is positive, which is probably a response to the 
ageostrophic u-component created as the front undergoes 
adjustment. At 48 hours the v-field has a small positive 

and may be inertial motion. 



B. SOMMAB Y -- CASE III 



The application of heating is able to generate and pre- 
serve a warm plume which overrides the cooler water in 
response to net Ekman transport. Some dcwnwelling is 

implied at the leading edge of the front. The isotherms 
show a "kirk ,, and an ageostrophic u-velocity is there. Seme 
of the surface-warmed water is mixed into the ambient pool. 
On the right "upstream" side, upwelling occurs as cooler 
water is brought up to replace that lost by transport. 
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Figure 34. Front 2 Case III 36- Hour U-Velocity. Cor.tcu: 
increment is 20 cm/s. 
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Figure 35. Front 2 Case III 24-Hour V-Velocity. Contour 
interval is 1G cm/s. 
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The mixed layer depth shallow s to a constant depth outside 
of the frontal area under the uniform heating, and the shal- 
lower peak existent within the front is advected with the 
front in a response to the mass transport. The across-front 
velocity generally shows a maximum at the surface within the 
front the direction of net Ekman transport, though reversals 
in direction occur with time as a result of ageostrophic 
re spcnses . 
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M ODE L RESULTS FOR FRONT 2 CASE IV 
A. ANALYSIS AND DISCUSSION 

The combined effects cf surface hearing and a negative 
wind stress act to push the upper 150 m of the front into 
the warmer water on the right side. Vertical mixing is pro- 
nounced as seen by the vertical orientation of the surfaced 
isotherms (Fig. 36) . This result is expected as convective 
mixing is coupled with wind stirring. Diffusion and hori- 
zontal mixing cause the frontal horizontal temperature gra- 
dient to decrease continually over the 48-hour integration. 
Here too, the entire field develops a downslope to the right 
as in the preceding three cases. 

The mixed layer is initially shallowed to around 40- m by 
the applied heating, and the shallow section located within 
the front breaks down and approaches the level of the over- 
all mixed layer (Fig. 37) . This result is expected based on 
previous cases which combined the action of convective mix- 
ing and wind stirring. 

The u-velccity field has an ageostrophic component at 
the left boundary as a response to the adjustment. 



101 



Depth (m) 



T at hour 36 




Fiqure 36. Front 2 Case IV 36-Hcur Temperature. Contour 
interval is 0.5 deg C. 
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Figure 37. Frcnt 2 Case IV 36-Hour Mixed Layer Depth. 
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Basically the u-vslocity decreases its maximum speed and 
weakens its horizontal gradient over the 48-hcur 
integration. The across-front v-velocity is in the direc- 
tion cf transport at 12 and 48 hours and extends from the 
surface, through the mixed layer depth, and down to a depth 
below 100 m. Once again, cross-f ror.tal mixing occurs. At 
24 hours, the v-field is reversed over the entire upper 2C0 
m. At 36 hours, the v-field exists only in the shallow sur- 
face layer, and, by 48 hours, the field returns to the posi- 
tive direction. Some type of oscillation is evident in the 
v field (refer to Figs. 38, 39, 40, and 41). At hour 12, 
the field is positive everywhere. At hour 24, it is neg- 
ative everywhere. At hour 36, it is negative in the shallow 
depths only. At 48 hours, it is again positive everywhere. 
This pattern occurred in the previous three cases, and it is 
strongly suggestive of seme type of oscillation in the 
v-component both in direction and in depth-extent. 

B. SUMMARY -- CASS IV 

The adjustment of the front under a negative wind stress 
and surface heating is affected by both convective mixing 
and wind-stirring in the vertical. Cross-f ror.tal mixing is 
evident in the horizontal. The front is "pushed back" into 



the warmer water on the right by the transport of coder 
water. The front diffuses over time. 

The mixed layer depth shallows outside of the front due to 
the surface heating, and the shallow portion within the 
front is eroded through nixing and approaches the depth of 
the far field. The u-velccity fellows an uneventful pattern 
of diffusion and a decrease in maximum speed as the front 
adjusts. The direction of the v-velocity correlates to the 
direction of net Ekman transport only at 12 and 48 hours. 
Between those times, an oscillatory pattern in direction and 
depth is observed. This is the same pattern which was 
observed in the v-fieids of Cases I, II, and III. The 
v^field indicates cross-f ront al mixing as it extends over a 
great depth. It is independent of the mixed layer depth. 
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F-qure 38. Front 2 Case IV 12-Hcur V-Velccity. Contour 
interval is 5 cm/s. 
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V at hour 24 




. Front 2 Case IV 24-Hcur V-Velocity. 
interval is 5 cm/s. 
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Figure 39 



Contour 



V at hour 36 




Fiaure 40. Front 2 Cass IV 36-Hcur V-Velccity. Contour 
interval is 5 cm/s. 
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Figure 41. Front 2 Case IV 48-Hcur V-Velccity. Contour 
interval is 5 cm/s. 
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XI. CONCLUSIONS AND RECOMM ENDATI ONS 

A. GENERAL 

This thesis has examined the unsteady response of two 
oceanic fronts to the local atmospheric forcings of wind 
stress and surface heat flux. The simplicity of Front 1 
enabled identification of basic processes stimulated by the 
forcings, processes which were not so clearly seen in Front 
2. Front 2 was a mere realistic simulation of an oceanic 
front, and, hopefully, it will serve as a stepping stone for 
further investigations intc fronts and testing of this fron- 
tal model. 

There are a few very general properties that exist for 
both Fronts 1 and 2, and which bear out our hypotheses 
stated in Chapter 2.B.: 

1) Wind direction plays a primary role in frontal 
adjustment. In the cases where warmer, less dense water is 
transported tc cccler, denser water, the structure of the 
front is preserved. The mixed layer depth is also pre- 
served, and little crcs s- frontal mixing occurs. 

2) In the cases where cooler, denser water is trans- 
ported to warmer, less dense water, the horizontal tempera- 
ture gradient is diffused, the mixed layer depth within the 
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front breaks dcwr. and deepens, and much cross-frontal fixing 
occurs. Ccnvective mixing occurs along with wind stirring. 

3) A uniform constant surface heat flux shallows the 
mixed layer outside of the frontal zone to a depth which is 
proportional tc the amount cf heat flux. The combination of 
a +1 dyne/cm^ wind stress and a surface heat flux of -0.004 
deg-cm/s mixes heat down to a depth of the initial mixed 
layer and leads tc stratification there. Stratification of 
the isotherms at the surface did not occur. For the cases 
in which the heat flux is applied, turbulent mixing is 
restrained in the vertical, which agrees with the results of 
Niiler (1S75) and Mellor and Durbin ( 1975) . 

4) In areas of the frontal zona where agsost rophic 
velocities are created, vertical mixing is increased. 

5) Mixed layers deepen rapidly within the first 12 
hours under wind stress alone as the fronts undergo initial 
adjustment. Eeycnd 12 hours, ad justments are confined tc an 
area around the the basic depth attained during the initial 
adjustment period. Niiler (1975), Mellor and Durbin (1975) 
and De Szceke (1980) all reported this phenomenon. 

6) A time-dependent pattern in the direction and 
depth-influence of the across-fronx velocity exists. It is 



probably due tc inertial oscillations created by the atmos- 
pheric forcing. 

It was surprising to see the downward bulge in the mixed 
layer at the cold-side boundary of both fronts in the wind 
stress only cases. It is unclear as to why this happens, 
though some speculative reasons ware provided. 

In Front 1, the formation of cores of positive v-veicc- 
ity which are symmetric about the mixed layer depth are 
intriguing. Inertial oscillations may provide an 

ex planaticn. 

The tilting down to the right of the entire temperature 
field in Front 2 in ail four cases is a major adjustment, 
and one which is independent of wind stress direction and 
heating. It would appear that this tilt is a response to 
the re-orientation of the oceanic current fields. 

B. K ECO M MENTATIONS 

There is both further investigation that can be done 
with the results of this thesis, as well as new work which 
is suggested by these results. The results of these eight 
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-dimensional version of this 
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he field would provide fut- 
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(1977) of the Maltese front 
point . 

nit (CPU) time was extensive 
They each averaged over four 



9 



hours of CFU tins to run, providing no stopping and restart- 
ing were required. Such a model is far from being useful 
operationally for real-time assessment of ocean structure to 
provide, for example, input for acoustic forecasting models. 
However, fcr basic research, this model appears to be a pow- 
erful tccl for deciphering the complexities of ocean frontal 
dynamics and thermodynamics. 
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